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Abstract 

Using  the  moving-singularity  finite  element  method  described  in  Part  I 
of  this  paper,  several  problems  of  dynamic  crack  propagation  in  finite  bodies 
have  been  analysed.  Discussions  of  the  effects  of  wave  interactions  on  the 
dynamic  stress-intensity  factors  are  presented.  The  obtained  numerical  results 
are  compared  with  the  corresponding  infinite  domain  solutions  and  other 
available  numerical  solutions  for  finite  domainsy| 

Introduction  / 

In  Part  I  of  the  present  paper  [l],  a  "moving  singular-element"  procedure 
has  been  presented  for  the  dynamic  analyses  of  problems  of  fast  crack  propagation 
in  arbitrarily  shaped  finite  bodies  with  linear  elastic  material  behaviour.  In 
this  procedure  a  singular-element,  within  which  a  large  number  of  analytical 
eigen-functions  corresponding  to  a  propagating  crack  are  used  as  basis  functions 
for  displacements,  may  be  translated  by  an  arbitrary  amount  A£  in  each  time 
increment  At  of  the  numerical  time-integration  procedure.  The  moving  singular- 
element,  within  which  the  crack-tip  has  always  a  fixed  location,  retains  its' 
shape  at  all  times,  while  the  mesh  of  "regular"  (isoparametric)  finite  elements, 
surrounding  the  moving  singular  element,  deforms  accordingly.  An  energy- 
consistent  variational  statement  was  developed,  as  a  basis  for  the  above  "moving 
singularity"  finite-element  method  of  dynamic  crack  propagation  analysis.  II 
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has  been  shown  [il]  that  the  present  procedure  leads  to  a  direct  evaluation 
of  the  dynamic  stress-intensity  factors. 

In  .the  present  Part  II  of  the  paper  several  numerical  studies  of  stationary 
as  well  as  propagating  cracks  in  finite  bodies,  are  presented.  These  studies, 
in  general,  fall  into  the  category  of  linear  elastic  dynamic  fracture  mechanics. 
These  studies,  in  addition  to  illustrating  the  efficiency  and  accuracy  of  the 
present  procedure,  also  shed  light  on  the  effects  of  stress-wave  interactions 
on  the  stress-intensity  factors  for  dynamically  propagating  cracks  in  finite 
bodies. 

The  presently  considered  examples  include:  (i)  static  extension  of  a 
central  crack  in  a  panel  from  a  non-zero  initial  length;  (ii)  self  similar, 
constant-velocity ,  propagation  from  a  finite  initial  length  of  a  central  crack 
in  a  finite  plane  body  subject  to  an  uniform,  time-independent,  tensile 
stresses  (normal  to  crack-axis)  at  the  edges  (which  problem,  is  analogous  to 
that  treated  by  Broberg  [2],  and  Rose  [3,4]);  (iii)  a  stationary  central  crack 
in  a  finite  plane  body  subject,  at  its  edges,  to  a  Mode  I  type  uniformly 
distributed  stress  with  a  Heaviside  step-function  time-dependence  (analogous 
to  the  problems  of  Baker  [5],  Sih,  Embley  and  Ravera  [6],  and  Thau  and  Lu  [7]); 
(iv)  a  problem  similar  to  that  in  (iii)  except  that  the  crack-tips  remain 
stationary  until  a  time  tQ,  where  upon  they  start  propagating  at  a  constant 
speed  (analogous  to  the  problems  studied  by  Freund  [8]);  and  (v)  constant- 
velocity  propagation  of  an  edge-crack  in  a  panel,  with  the  direction  of 
propagation  being  pare! lei  to  the  panel-edges  on  which  uniform  displacements, 
normal  to  the  edges,  are  prescribed  (analogous  to  the  problem  treated  by 
Nilsson  (')]).  All  the  references  [ 2-8 j  deal  with  unbounded  bodies,  except 
(9)  which  deals  with  a  finite  height,  but  infinite  width  strip.  Thus  the 
presently  obtained  results  for  finite  bodies  are  compared  with  those  in  [2-i] 
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and  the  effects  of  finiteness  of  the  domains  are  discussed.  Comparisons  of 
the  present  results  with  the  numerical  results  of  other  investigators,  where 
available,  are  also  presented  and  discussed. 

In  the  following  we  present  results  for  each  of  the  problems  cited  above. 

(i)  Static  Crack  Extension 

To  test  the  accuracy  of  the  present  method  of  "moving  singular-elements", 
first  a  static  problem  of  a  central-cracked  square  panel  [2L(length)  =  2W(width)J, 
subjected  to  uniform  tension  at  edges  parrellel  to  the  crack  axis,  was  solved 
to  obtain  the  static  stress-intensity  factor  as  a  function  of  the  current  crack 
length,  Z.  Thus,  in  the  finite  element  development  given  in  Part  I  of  this 
paper  (1  ],  velocity  and  acceleration  effects  were  ignored.  The  eigen-functions 
embedded  in  the  singular-element  reduce,  when  v  =  0,  to  the  well-known 
Williams'  eigen-functions  as  shown  in  Appendix  A  of  [l].  Starting  from  an 
initial  crack  length  value  of  ZQ  =  0.2W,  the  singularity  element  was  successively 
translated  (in  Mode  I  growth  sense),  statically  (with  v  =  0)  in  increments  of 
AE  =  0.005W  until  a  final  value  of  crack  length  Z  =  0.5W  is  reached.  During 
the  above  series  of  calculations,  the  externally  applied  uniform  tension  was 
held  constant.  At  each  current  crack-length  level,  the  stress-intensity 
factor  is  computed  directly  as  an  unknown  from  the  finite  element  equations,  as 
described  in  [1  ].  From  the  normalized  stress-intensity  factor  solutions  shown 
in  Fig.  1,  it  is  seen  that  the  present  results  agree  excellently  with  those 
reported  by  Isida  Do]-  The  normalized  results  shown  in  Fig.  1  may  be  viewed 
as  correction  factors  for  static  stress-intensities,  due  to  the  finite  size  of 
the  oanel . 

,,e  note  chat  the  symbols  (  t)  in  Fig.  1,  as  well  as  in  all  the  subsequent 
Figures,  denote  the  current  crack  length  (Z/W)  where  the  regular  elements 
surrounding  the  moving  crack-element  were  readjusted  as  described  in  Fig.  ? 
of  11  1. 
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The  problem  is  that  of  a  centrally  cracked  square  panel  (L  =  W  =  40mm) 
with  properties:  y  {shear  modulus)  =  2.94  x  10  N/rmr;  v(  Poisson's  ratio) 

=  0.286;  p  (  Mass  density)  =  2.45  x  10^  K^m^.  A  time-independent  tensile 
stress  was  assumed  to  be  acting  at  the  edges  of  the  specimen  parrel  lei  to  the 
crack-axis.  The  crack  is  assumed  to  open  from  an  initial  length  (£Q/W)  =  0.2 
and  to  grow  symmetrically  with  a  constant  velocity,  v.  This  problem  may  be 
considered  to  be  similar  to  that  treated  by  Broberg  [2]  except  that  Broberg 
treated  an  infinite  body  with  a  crack  which  opens  from  a  zero  initial  length. 
The  problem  was  analyzed  for  four  different  values  of  v,  namely,  (v/Cs)  =  0.2, 
0.4,  0.6  and  0.8  respectively,  where  the  shear  wave  speed  for  the  present 
problem  is  C$  =  3.4641  x  10^  mm/sec.  The  dilatational  and  surface  (Rayleigh) 
wave  speeds,  Cd  and  CR,  respectively,  are  such  that:  (Cd/Cs)  =  1.8266,  and 
(CR/C$)  =  0.9238. 

In  all  the  four  considered  cases  of  (v/Cs)  ratio,  the  increment  of  crack 

growth  in  each  step,  AZ,  was  kept  the  same,  at  the  value:  (AZ/W)  =  0.005. 

Thus,  in  each  of  the  considered  (v/Cs)  cases,  the  time  integration  step.  At, 

changes  according  as:  [(v.At)/z]  =  0.005.  The  finite  element  mesh  used,  at 

the  initial  crack  length  in  each  of  the  four  cases,  is  shown  in  the  inset  of 

Fig.  3.  In  this  Figure,  as  well  as  in  Figs.  (1,6  and  9),  the  singular-element 

near  the  crack-tip  is  identified  by  hatched  markings. 

As  noted  in  the  review  article  by  Rose  [3],  the  dynamic  stress  intensity 

factor  K  may  be  expressed  as  the  product  of  a  velocity  factor  k(v)  and  a 
★ 

static  factor  K  ;  thus, 

K  =  k(v)  K*  (1) 

The  "static  factor"  K*  depends  on  the  current  length  of  the  crack,  the  applied 
load,  the  history  of  crack  extension,  but  not  on  the  instantaneous  crack  speed 


As  also  discussed  in  [3],  K  is,  in  general,  not  equal  to  the  static  stress- 
intensity  factor,  Ks,  for  a  stationary  crack  of  the  same  length  as  the 
moving  crack.  The  analytical  expression  for  !<*°°  (in  an  infinite  body  sub.i'ic !; 
to  uniform  stress  normal  to  the  crack  axis),  as  a  function  of  the  current 
crack  length  is  given  by  Eshelby  [11]  and  in  [3],  as: 

K*~  -  (2/tt)55  f  a  (x)/(E  -  *)h  dx  (2) 


where,  a  (x)  =  a  |x|/(x2  -  E2)"5  ;  |x|  >  E 


In  Eqs.  (2  and  3),  a  is  the  applied  tensile  stress  at  infinity,  2Eq  is  the 
initial  crack  length,  and  x,y  are  cartesian  coordinates  centered  such  that  x 
=  +  E(t)  denote  the  current  crack-tips  and  y  is  normal  to  crack-axis.  Thus 
°yy  (x)  in  Eq.  (2)  is  the  initial  distribution  of  stress  along  the  axis  of 
the  crack,  prior  to  its'  propagation.  Eshelby's  results  for  the  integral  in 
Eq.  (2)  is: 

K*°°  =  [o(nEo)1*  ( 2/ it )  {(2+0  E-F}]/[1  +  U/2)]5*  (0 

where  £  =  (E  -  Eq)/Eq-,  and  E  and  F  are  the  complete  elliptic  integrals  of  the 

I- 

2nd  and  1st  kind;,  respectively,  with  the  modulus  [£/(2+£)]  4.  On  the  other 
hand,  the  static-stress  intensity  factor  K~  for  a  crack  of  length  2E  in  an 


infinite  domain  is, 


tC  =  o(  ttEP 


Thus,  in  general,  K  <  k”.  It  is  also  noted  that  the  Eq.  (4)  for  K  is 
valid  only  until  the  time  that  disturbances  from  one  crack-tip  reach  the  oth'-r 
(moving)  crack-tip.  However,  for  the  case  of  a  crack  growing  sel f-similarly 
from  a  zero  initial  length  (r  =  0)  with  a  constant  velocity  v,  as  in  the 
problem  studied  by  3roberg  [21,  disturbances  from  one  crack  tip  influence  at 
all  times  the  other  moving  crack  tip  if  2v  <  ch  ,  and  K*10  t  £  at  all  times. 


The  normalized  dynamic  stress-intensity  factor  solution  for  the  present 
problem  of  a  crack,  in  a  finite  square  panel,  propagating  self-similarly  at 
a  constant  velocity  (v/C$)  =  0.2, is  shown  in  Fig.  2.  Note  that  the  normaliza¬ 
tion  is  such  that  a  value  of  unity  represents  the  normalized  stress  intensity 
factor  K“  for  a  crack  of  length  2Z(Z  =  I  +  vt)  in  an  infinite  solid  with  time- 
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independent  remote  tension  a.  Also  shown  in  Fig.  (2)  are:  (i)  the  finite-size 
correction  factor  F  =  K^/ctCttE)^  for  the  static  stress-intensity  factor  in 
the  present  finite  domain;  (ii)  the  normalized  "static"  factor  K  °°/a(iT>:)'2  , 
as  calculated  from  Eq.  (4),  and  (iii)  the  velocity  factor  k(v)  for  ( v/Cs )  =  0.2, 
as  given  by  Broberg  D3  .  The  effect  of  the  finiteness  of  the  domain  on  K 
may  be  accounted  for,  approximately,  by  including  a  finite-size  correction  factor 
in  the  initial  stress-distribution  at  ZQ,  in  Eq.  (3)  above.  Thus  it  appears 
that  one  may  write,  approximately ,  that  K*f  =  F( E0) K  °°.  For  the  present 
case  of  ( E0/W)  =  0.2,  the  finite  correction  factor  F( )  =  1.055,  as  seen  from 
Fig.  1. 

Also  marked  in  Fig.  2  are  several  specific  instants  of  time  (or  equivalently), 
the  corresponding  values  of  S)  with  the  notations:  (a)  Dc,  Sc  and  Rc  are, 
respectively,  the  times  taken  by  the  dilatational ,  shear  and  Rayleigh  waves 
to  travese  one  crack-width,  (b)  DCDC  (or  RCRC)  denotes  Dc  (or  Rc)  plus  the 
time  for  the  first  rescattered  dilatational  (or  Rayleigh)  waves  to  travel  one 
crack-width;  (c)  (or  SCS^)  denotes  the  time  taken  for  the  dilatational 
(or  shear)  waves  emanated  by  one  crack-tip  to  be  reflected  by  the  nearest 
free-boundary  and  reinteract  with  the  crack-tip  in  question.  The  above  times 
are  calculated  from  the  continuum  relations  for  the  respective  wave  speeds. 

However,  it  should  be  borne  in  mind  that  a  "consistent-mass"  representation 
is  used  in  the  present  finite  element  method. 


For  the  case  of  (v/C$)  =  0.2,  it  is  seen  from  Fig.  2  that  the  computed 
normalized  dynamic  stress-intensity  factor  correlates  excellently  with  the 
values  given  by  F(Eo)K*“k( v)/a(7rE)'s  until  roughly  the  time  denoted  by  Rc> 

At  longer  times,  ie.,  at  the  times  greater  than  R  and  D  Df,  the  computed 
normalized  dynamic  Kj(t)  appears  to  correlate  excellently  with  the  values 
given  by Jb  K^k(v)  /  a(iTE)'5]  where  0  is  a  constant.  It  is  interesting  to 
observe  that,  for  the  present  problem,  this  constant  0  appears  to  be  equal 
to  g(£rc)/F(Erc)  where  F(ERc)  is  the  finite  correction  factor  in  the  static 
stress-intensity  for  a  crack  df  length  equal  to  the  current  length,  ERc,  in 
a  dynamic  problem,  at  which  the  event  Rc  (as  defined  earlier),  occurs;  where¬ 
as,  G(W  1S  likewise,  the  ratio  [f(Eq)K  “/o/tteJ  at  E  =  ERc.  The  event  R£ 
is  seen  to  occur  at  the  time,  t  ==  (2EQ)/(CR-v) ,  and  thus,  ERc  =-E  +  vt. 

Both  t  and  (eRc/W)  decrease  as  EQ  decreases,  for  given  CR  and  v.  Thus,  for 
smaller  values  of  EQ  (and/or  smaller  values  of  v),  the  ratios  F(ERc),  G( 
and  hence  0  tend  to  a  value  of  unity.  Thus,  for  cracks  propagating  from 
initial  lengths  such  that  (Eq/W)  «  1,  it  appears  that  at  the  times  greater 

than  R£  and  DcD^  and/or  after  the  crack  has  grown  dynamically  to  a  few  times 
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its*  initial  length,  the  static  factor  K  approaches  the  static-stress- 
in tensity  factor  for  the  current  crack  length  in  the  finite  body,  namely,  !<  . 

In  connection  with  the  presently  computed  resul ts  shown  in  Fig.  2,  it 
should  be  noted  that  the  crack  velocity,  v,  was  taker,  to  be  zero  at  the  initial 
crack  length  E  .  It  is  assumed  that  the  crack-tip  accelerates  to  a  velocity, 
(v/Cs)  =  0.2,  during  the  first  time  increment.  At.  The  convergence  of  the 
present  numerical  results  to  the  analytically  predicted  ones,  at  small  times, 
t  <<  U  in  the  present  case  of  (v/C  )  =  0.2,  could  have  been  studied  by  altar  in 
this  time  step  At  in  which  the  crack-tip  accelerates  from  (v/Cj  =  0  to  (v/r  1 
=  0.2.  However,  this  was  not  attempted. 


The  computed  normalized  dynamic  stress-intensity  factor  solution  for  the 
case  (v/C$)  =  0.4  is  shown  in  Fig.  3,  wherein  the  times  Dc,  Sc,  DcDf  and  R  , 
as  defined  earlier,  are  also  marked.  Once  again,  it  is  seen  that  until  signifi¬ 
cant  interaction  of  the  waves  from  the  other  crack-tip  and  the  free  surface 
takes  place  (ie.,  for  time  t  <  Rc  or  0CD^),  the  computed  dynamic  K-factor  for 
the  finite  body  correlates  excellently  with  the  value  predicted  by  the  approximate 
function:  ,  £K*fk(v)j. 

Finally,  the  results  for  the  cases  (v/Cs)  -  0.6  and  0.8,  respectively, 
are  shown  in  Fig.  4,  wherein  only  the  times  Dc,  and  DcD^  are  also  marked. 

The  values  of  R  and  S  are  greater  than  the  time  for  which  the  solution  is 
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obtained.  Once  again,  it  is  seen  that  the  computed  dynamic  K-factor  correlates 
well  with  the  approximate  prediction,  K*^k(v)  .  However,  the  convergence  cv 
the  computed  solution,  to  that  analytically  predicted,  is  slow,  at  these 
higher  crack-speeds.  A  possible  reason  for  this  may  be  the  initial  conditions 
at  EQ  used  in  the  present  study,  as  explained  eariler. 

The  crack-mouth  opening  displacements  at  various  instants  of  time  (or 
equivalently,  at  corresponding  crack-lengths),  for  the  case  of  (v/Cs)  =  0./l , 
are  shown  in  Fig.  5.  Also  shown  in  Fig.  5  are  the  corresponding  analytical 
results  by  Broberg  [2],  who  consideres  constant  velocity  crack-propagation 
starting  from  a  zero  initial  crack-length.  An  excellent  correlation  between 
the  present  results  and  those  of  [2]  is  noted.  A  similar  correlation  was 
also  noted  for  the  other  considered  cases  of  (v/C$)  ratio,  but  are  not  shown 
here. 

( i i i )  Stress-Wave  Loading  of  a  Stationary  Crack  in  a  Finite  Body 

The  problem  is  that  of  a  rectangular  panel  [(w/L)  =  2.6]  with  a  centrally 
located  crack  of  length  (Hq/W)  =  (3/13).  The  material  properties  are  taken 
to  be:  p(shear  modulus)  =  2.94  x  10^  N/nm^ ;  v  =  0.286;  and  p  =  2.45  x  10J 
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kg/m  .  Uniformly  distributed  uniaxial  tensile  stresses,  with  a  Heaviside  step- 
function  time-dependence,  were  assumed  to  act  at  the  edges  of  the  panel  parrel  lei 
to  the  crack-axis.  The  crack  is  assumed  to  be  stationary  under  the  action  of 
this  time-dependent  loading. 

Due  to  symmetry,  only  a  quadrant  of  the  panel  is  modeled  by  finite  elements, 

as  shown  in  Fig.  6.  Also  marked  in  Fig.  6  are  specific  instants  of  time, 

calculated  by  using  continuum  wave  speeds,  with  the  notations:  (a)  Dm  is 

the  time  taken  by  the  dilatational  waves  to  travel  the  distance  from  the 

boundary,  where  time-dependent  tractions  are  applied,  to  the  crack-tip;  (b) 

DmDc’  DmSc  and  DmRc  are>  respectively,  equal  to  D  plus  the  time  taken  by 

the  first  scattered  dilatational,  shear,  and  Raleigh  waves  to  reach  from  one 

crack-tip  to  the  other;  (c)  DmDcDc  is  equal  to  DmDc  plus  the  time  taken  by 

the  first  rescattered  dilatational  wave  to  travel  one  crack-width;  (d)  D  0  D. 

'  m  c  v 

is  equal  to  Dm  plus  the  tirte  taken  by  the  scatter 2.'.  longtitudina!  -paves  to 
travel  from  the  crack-tip  to  the  nearest  free  boundary  surface  and  back  to 
the  same  crack-tip;  (e)  DmDn(  is  the  time  taken  by  the  dilatational  waves  to 
travel  the  length  of  the  bar,  reflect  from  tiie  boundary  surface  on  the 
opposite  and  return  back  to  the  crack-tip;  and  (f)  DmDn)Sc  and  DmDnR 

are,  respectively,  equal  to  DmDm  plus  the  time  for  the  rescattered  dilatation?!  , 
shear  and  surface  waves  to  travel  one  crack-width . 

The  presently  computed  normalized  dynamic  stress-intensity  factor  solution 

is  shown  in  Fig.  6.  Also  sho wn  in  Fig.  6  are  the  analytical  solutions  by 

Baker  [5J ,  Sih,  Embley  and  Ravera  [6],  for  infinite  domains,  and  the  numerical 

solution  by  Aoki  et  al  [rij  for  a  finite  domain  identical  to  the  one  consider?-' 

here.  During  the  time  interval  D  to  D  D  ,  when  no  wave  interaction  takes 

in  lil  c 

place,  the  results  for  the  present  problem  must  agree  with  the  results  of 
Baker  [s]and  this  can  be  seen  to  be  the  case  with  the  present  results.  The 
solution  at  longer  times  is  found  to  be  in  good  agreement  with  that  of  Sih  et  al 


[6].  It  is  noted  that  the  overshoot  in  the  k-factor  at  the  time  instant  DmRc 
in  the  present  solution,  as  compared  to  the  solution  by  Sih  et  al  [6],  is 
analogous  to  that  in  a  recent  solution  by  Kim  03] .  It  is  also  seen  from  Fig. 
6,  that  the  present  solution  is  higher  than  that  of  Aoki  et  al  [12]  at  all 
times.  However,  the  solution  by  Aoki  et  al  appears  to  be  lower  than  that  by 
Baker  [5]  even  for  times  less  than  DmDc. 

Finally,  it  is  of  interest  to  note  that  even  though  the  time  for  the 

longitudinal  wave  to  arrive  from  the  loaded  boundary  to  the  crack,  as  computed 

from  continuum  wave  speeds,  is  Dm  as  marked  in  Fig.  6,  a  nonzero  stress- 

intensity  is  observed  at  the  crack-tip  even  at  times  lower  than  D  in  the 

m 

finite  element  solution.  This  is  due  to  the  inertia-coupling  that  exists 
between  the  finite  element  nodes  (especially  those  of  the  singular-element) 
when  a  consistent  mass-matrix,  as  in  the  present,  is  used. 

(iv)  Crack-Propagation  at  Constant  Speed:  Stress -Wave  Loading 

We  consider  the  problem  wherein  the  geometry,  material  properties  and 
the  time-dependent  loading  are  all  identical  to  those  described  under  the  case 
(iii)  above.  In  this  problem,  the  crack  with  an  initial  length  of  Uq/W) 

=  3/13  remains  stationary  until  a  time  t  =  4.4  ps  and  then  propagates  with 
a  constant  velocity  v  =  1000  m/sec.  The  finite  element  breakdown  at  the 
initial  crack  length,  is  identical  to  that  in  Fig.  6.  In  modeling  the 
crack  propagation,  the  regular  elements  are  periodically  readjusted  as 
indicated  in  Fig.  2  of  [l].  The  instants  of  time  (or  equivalently  the  value 
of  7.)  at  which  this  readjustments  are  done,  are  marked,  by  (t)  in  Fig.  7. 

Also  marked  in  Fig.  7  are  several  specific  instants  of  time,  calculated 
by  using  continuum  wave  speeds,  with  the  notations:  (a)  D  is  the  time  for 
the  dilatational  waves  to  travel  from  the  boundary,  where  a  Heaviside  step 
function  tension  is  applied,  to  the  stationary  crack-tip.  (b)  DmDc,  DmSc 
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and  D  R  are,  respectively,  equal  to  D  plus  the  time  taken  by  the  di latational , 

shear  and  surface  waves  first  scattered  by  the  stationary  crack-tip  to  reach 

the  other  crack-tip;  (c)  Dc,  Sc,  Rc  are,  respectively,  equal  to  tQ  (when  the 

crack-tips  begin  to  propagate)  plus  the  time  taken  by  the  dilatational ,  shear 

and  surface  waves  emanated  by  crack-tip  at  tQ  to  reach  the  other  crack-tip; 

(d)  D  D  is  the  time  taken  by  the  longitudinal  waves  to  travel  the  length 
''mm 

of  the  panel  to  the  opposite  side  and  return  to  a  (moving)  crack-tip;  (e) 

D  D  is  the  time  taken  for  the  dilatational  waves  emanating  from  a  crack-tip 
cm 

at  tQ  travel  to  the  nearest  boundary  (which  in  this  case  is  the  one  where 
tractions  are  applied)  and  back  to  the  same  crack-tip,  (f)  is  equal 

to  Dm  plus  the  time  taken  by  the  dilatational  waves  first  scattered  by  the 
stationary  crack-tip  to  travel  to  the  nearest  free-boundary  and  back  to  the 
same  crack-tip,  which  is  now  propagating;  (g)  DcDf  is  equal  to  tQ  plus  the 
time  taken  by  the  dilatational  waves  emanating  from  the  crack-tip  at  tQ  to 
travel  to  the  nearest  free  boundary  and  back  to  the  same  crack-tip. 

The  presently  computed  results  for  the  dependence  of  the  dynamic  stress 
intensity  factor  on  time  are  indicated  in  Fig.  7,  along  with  the  comparison 
(analytical)  results  by  Freund  [8]  and  Baker  [5],  and  the  (numerical)  results 
by  Aoki ,  et  al  [12]  .  As  noted  earlier,  the  correlation  of  the  present  stationary 
crack  results  with  those  of  Baker  [5],  until  wave  interaction  takes  place,  is 
excellent.  In  the  case  when  the  crack  propagates  suddenly  at  the  time  t 
(=  4.4  us  in  Fig.  7),  the  infinite  domain  results  by  Freund  [3]  are  seen  to 
correlate  excellently  with  the  present  results,  until  significant  interaction 
of  waves  emanating  from  one  crack-tip  with  the  other  crack-tip  takes  place. 

It  is  noted  that  in  [d]  the  crack-tip  velocity  is  assumed  to  change  from  zero 
to  "v"  in  zero  time,  whereas,  in  the  present  numerical  study,  this  transition 
is  assumed  to  take  place  over  a  finite  time-step  At.  For  a  closer  comparison 


with  the  results  of  [8],  it  would  be  interesting  to  vary  the  size  of  At  over 
which  the  above  transition  occurs;  but  this  is  not  pursued  in  the  present 
study.  .It  is  seen  that  the  comparison  results  by  Aoki  et  al  [12]  are  somewhat 
lower  than  those  in  [8]  even  until  the  time  that  the  solution  in  [8]  may  be 
considered  valid.  Further,  the  results  in  [lij  are  higher  than  the  present, 
after  this  time. 

The  crack-face  opening  displacements  at  various  time  intervals  for  the 
propagating  crack,  as  well  as  similar  results  at  corresponding  times  for  a 
stationary  crack  are  shown  in  Fig.  8.  It  is  interesting  to  note  that  at 
sufficient  distances  away  from  the  propagating  crack-tip,  the  crack-mouth 
opening  displacements  are  nearly  the  same  for  the  stationary  as  well  as. 
propagating  cracks. 

(v)  Constant  Velocity  Crack-Propagation  in  a  Strip  with  Prescribed  Boundary- 
Displacements 

The  problem  considered  is  that  of  the  constant  velocity  propagation  of  an 
edge  crack  in  a  square  sheet  whose  edges  parrellel  to  the  direction  of  crack- 
propagation  are  subject  to  uniform  displacements  in  the  direction  normal 
to  that  of  crack-propagation.  This  problem  is  analogous  to  that  treated  by 
Nilsson  [9]  who  obtained  an  analytical  solution  for  the  steady-state  stress- 
intensity  factor  for  the  constant  velocity  propagation  of  a  semi-infini te 
crack  in  a  finite-height  (normal  to  crack-axis),  infinite-width,  strip. 

In  the  present  problem  the  following  geometry  and  material  parameters 

are  used:  (2h/U)  =  1.0  (see  inset  of  Fig.  9),  v  =  0.286;  u  =  2.94  x  1010 
2  3 

and  0  =  2.45  x  10  Kg/m  .  The  crack  is  assumed  to  start  to  propagate  from 
an  initial  length  of  (Z  /h)  =  0.4. 

Three  different  cases  of  constant  velocity  propagation,  (v/C$)  =  0.2; 
0.40,  and  0.60,  respectively,  are  considered. 
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The  results  for  the  dynamic  stress-intensity  factor  for  the  case  (v/Cs) 

=  0.2  are  shown  in  Fig.  9  along  with  the  analytical  solutions  by  Nilsson  [9J. 

The  results  in  Fig.  9  are  normalized  with  respect  to  the  plane-strain,  static 
(v  =  0)  stress  intensity  factor  for  the  semi-infinite  crack  in  a  finite-heiaht, 

1  O 

infinite-width  strip,  namely,  K~  =  u2E  /  h^(lV)  .  Also  shown  in  Fig.  9 

are  certain  specific  instants  of  time,  with  the  notations:  (a)  DcD^-j , 

ScSfi ,  are,  respectively,  the  times  taken  by  the  dilational,  and  shear  waves, 

emanated  by  the  moving  crack-tip  to  be  reflected  from  the  nearest  free-boundary 
and  travel  back  to  the  crack-tip;  (b)  DCD^2  is  the  time  taken  by  the  dilatational 
waves  emanated  from  the  crack-tip  to  be  reflected  by  the  second  free-boundary 
(x  =  W)  and  travel  back  to  the  crack-tip. 

It  is  seen  from  Fig.  9  that,  the  correlation  between  the  present  and 
Nilsson's  [9]  results  is  excellent.  It  may  be  of  interest  to  note  that  in 
an  analysis  using  the  "node-release"  technique,  Mai  luck  and  King  ju]  found, 
for  the  case  of  (v/Cs)  =  0.325,  that  their  computed  results  for  the  energy- 
release  rate  for  a  similar  problem  were  about  20 %  lower  than  that  predicted 
by  Nilsson  [9],  as  steady-state  conditions  are  reached.  This  may  indicate 
the  relative  efficiency  of  the  present  numerical  method  as  compared  to  the 
familiar  "node-release"  techniques. 

Finally  the  computed  results  for  the  cases  ( v/C s )  =  0.4  and  0.6  are  shown 
in  Figs.  10  and  11  respecti vely.  The  times  noted  in  Figs.  10  and  11  have  the 
same  meanings  on  those  indicated  in  Fig.  9.  From  Figs.  10  and  11  it  is  once 
again  seen  that  the  present  results  agree  excellently  with  those  in  ^9j. 

Closure 

The  procedure  of  a  "moving  singular  element",  within  which  a  large 
number  of  eigenfunctions  for  a  propagating  crack  are  embedded,  has  been 
applied  to  study  several  problems  of  dynamic  crack  propagation  in  finite 


bodies.  The  numerical  results  have  been  found  to  correlate  well  with  the 
available  analytical  solutions,  for  corresponding  problems  in  infinite  do¬ 
mains,  during  the  time  for  which  these  analytical  solutions  may  be  considered 
as  valid.  The  computed  solutions  beyond  these  times,  and  the  knowledge  of  the 
times  involved  for  wave-interaction  in  finite  bodies,  indicate  both 
qualitatively  and  quantitatively  the  effects  of  stress-wave  interactions  on 
dynamic  stress-intensity  factors  for  cracks  propagating  in  finite  bodies. 

The  use  of  the  presented  numerical  procedure  in  the  simulation  of  experi¬ 
mental  data  from  dynamic  fracture  test  specimens,  as  well  as  in  predicting 
crack-propagation  history  in  dynamically  loaded  cracked  bodies,  are  the  sub¬ 
jects  of  a  forthcoming  paper. 
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Figure  Captions 


Fig.  1  Calculation  of  Static  Kj  factors  in  a  center-cracked  tension  specimen 

by  the  present  "moving-singularity"  method:  AE  =  0.005.  (+) 

indicates  the  current  crack  length  when  the  regular  element  are 
readjusted  as  shown  in  Fig.  2  of  [l]. 

Fig.  2  Normalized  Dynamic  Stress  intensity  factor  for  a  crack  starting  from 
a  finite  initial  length  and  propagating  with  constant  velocity, 
v/Cs  =  0.2. 

Fig.  3  Normalized  Dynamic  Stress  intnesity  factor,  ( v/C$ )  =  0.4. 

Fig.  4  Normalized  Dynamic  Stress  intensity  factors  for  the  cases  (v/C  )  = 

0.6  and  0.8. 

Fig.  5  Crack-face  Displacements  at  various  time-instants,  for  a  crack 
propagating  at  constant  velocity,  ( v/Cs )  =  0.4. 

Fig.  6  Time-Dependence  of  dynamic  stress-intensity  factor  for  a  center-cracked 
rectangular  plate  subjected  to  Heaviside  step-function  normal  stress. 

Fig.  7  Time-Dependence  of  dynamic  stress-intensity  factor  for  a  center-cracked 
rectangular  plate  subject  to  a  step-function  normal  stress.  Crack 
remains  stationary  until  tQ  and  then  propagates  with  a  constant  velocity. 

Fig.  8  Crack-face  profiles  for  a  stationary  as  well  as  propagating  crack  at 
various  times. 

Fig.  9  Dynamic  stress-intensity  factor  for  an  edge  cracked  square  sheet 
subject  to  constant  normal  displacement:  (v/C$)  =  0.2. 

Fig.  10  Dynamic  stress-intensity  factor  for  an  edge  cracked  square  sheet 
subject  to  constant  normal  displacement:  ( v/Cs )  =  0.4. 

Fig.  11  Dynamic  stress-intensity  factor  for  an  edge  cracked  square  sheet 
subject  to  constant  normal  displacement:  (v/Cs)  =  0.6. 
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